Probing the magnetic field with molecular ion spectra. 
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Abstract 

Observations of the effect of the magnetic field on its environment are usually 
achieved with techniques which rely on the interaction with the spin of the particles 
under study. Because of the relative weakness of this effect, extraction of the field char- 
acteristics proves to be a most challenging task. We take a totally different approach 
to the problem and show that the manifestation of the magnetic field can be directly 
observed by means of a comparison of the spectra of molecular ions with those of neu- 
tral molecules. This takes advantage of the strong cyclotron interaction between the 
ions and the field, but requires the presence of flows or turbulent motion in the gas. 
We compare our theory to data obtained on the OMCl, 0MC2, 0MC3 and DR210H 
molecular clouds. 

Subject headings: ISM: cloud — ISM: magnetic field — ISM: molecules 

1. Introduction. 

The suspected effects of the magnetic field in the process of star formation are well documented 
in the literature. They may occur during the initial collapse of clouds, by providing frictional 
support through collisions between ions and neutrals, or through the establishment of a pressure 
associated with long wavelength magnetohydrodynamic waves arising from small scale fluctuations 
in the field, or again as source of magnetohydrodynamic forces at the center of the molecular outflow 
phenomenon (Shu et al. 1987; Mouschovias 1991a,b; Crutcher et al. 1993). The need for accurate 
measurements of the magnetic field and its effects on the environment cannot be understated. This 
is, however, a formidable task. Up to now, the observed manifestations of the magnetic field, in star 
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forming regions, consist of Zeeman line broadening measurements (which are best carried out at 
low frequencies, due to the increasing dominance of the Doppler width over Zeeman splitting with 
increasing frequency) and in the magnetic alignment of the dust grains known to exist in molecular 
clouds. 

Despite the weakness of the interaction between the field and the observed molecular species 
and the numerous difficulties inherent to the technique, Zeeman measurements provide us with the 
only way of measuring the intensity of the magnetic field (more specifically its component normal 
to the plane of the sky (Crutcher et al. 1993, 1999)). On the other hand, a measure of the degree of 
alignment of the grains can be obtained from the level of linear polarization detected from continuum 
emission emanating from molecular clouds. The polarization is believed to be caused by anisotropic 
emission from dust grains. Recent studies (Hildebrand et al. 1999; Draine & Weingartner 1996) 
reveal that the intensity of the magnetic field has very little to do with the level of polarization 
induced by the dust. The field is needed to align the grains and the level of polarization is basically 
a function of the spin imparted to the individual grains by some agent, most likely the radiation 
field. The thermally generated spin is known to be too small to explain the observed effect (Draine 
& Weingartner 1996). 

In this paper, we present theoretical and observational evidence for a different effect which 
demonstrates the presence of the magnetic field. Its signature resides in the line profiles of molecular 
ion species at millimeter and submiUimeter wavelengths. It will be shown that, under the right 
conditions, even a weak field (~ 10 jiG) will cause ions to exhibit narrower profiles and a suppression 
of high velocity wings when compared to the lines of neutral species. 

There are of course several possible effects that could confuse the issue, including optical depth 
differences between species and lack of spatial coexistence. We deal with such problems in part in 
this paper and in more depth in a subsequent one. 

2. Ion spectra versus neutral spectra. 

Probably the first thing we should point out is that the lines observed in molecular clouds are 
usually many times broader than their thermal width. It is generally assumed that this is due to the 
presence of turbulence in the interstellar medium (Zuckerman & Evans 1974; Falgarone & Phillips 
1990). This will be the point of view adopted here. Turbulence is characterized by vortex streaming, 
the motions taking place in eddies of different sizes (Tennekes & Lumley 1972). In what follows, we 
will concentrate on a small portion of a given eddy so that the region considered is small enough 
that motions of the local fiow can be approximated as being linear, with the understanding that on 
a larger scale the motion exhibits vorticity. This analysis can be applied without approximation to 
cases involving linear fiows (e.g., jets, outfiows, . . . ). 

The effect to be discussed is the tendency for ions to be forced into gyromagnetic motion about 
the magnetic field direction, rather than following the streaming fiow of the general neutral motion. 
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It is necessary to discuss the ion-neutral collisions and we start with a simplified approach. We 

will assume that we are dealing with a weakly ionized plasma, where the neutral flow is mainly 
composed of molecular hydrogen; we therefore assume a molecular mass number of = 2.3 for its 
constituents. 



2.1. Approximate solution - head on collisions. 

The equation of motion of a single ion subjected to a flow of neutral particles in the presence 
of a magnetic field is given by: 



^ = — VXB + Fe, (1) 
at rriiC 

where mi, v, B are the mass of the ion, its velocity and the magnetic field. Fc is the force (per 
unit mass) on the ion due to the collisions with neutrals. The nature of the interaction during a 
collision can be quite complicated and how the momentum is transfered between the particles needs 
to be handled carefully. But in order to get an idea of the behavior of the ion we will initially assume 
that all the collisions are perfectly elastic and "head on". Accordingly, the force of interaction can 
be approximated by: 



Fe = -2^(V-V") Ht-Tm) 

rrii ^-^ 

m=— oo 



(2) 



where /i, v", r„j are the reduced mass, the neutral flow velocity and the different times at 
which collisions randomly occur. It will make things easier if we break up the velocities into two 
components, one parallel to the mean magnetic field (vy and vjj') and the other perpendicular (v_l 
and v"). In steady state conditions, it can be shown that the mean and variance of the velocity 
components are given by (see the appendix for details): 
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(iOr) is what we will call, for reasons to be discussed later, the relaxation rate. {uJ^), fc are the 
mean ion gyrofrequency vector and the (mean) collision rate. 

Equation (4) represents the drift that the ion can possess in relation to the neutral flow and it 
is at the heart of the effect that we are now studying. Equation (6) gives a measure of the gyration 
amplitude of the ion around a guiding center in the region occupied by the magnetic field. We have 
also allowed for an inherent velocity dispersion in the neutral flow as can be seen by the presence 



From equations (3)-(7) it can be deduced that, as could be expected, the ion completely follows 
the flow when the latter is aligned with the field (v^ = 0). More interesting, however, is the ion 
behavior when the neutral flow is perpendicular to the field (v||^ = 0). In such cases, the following 
observations can be made: 

• for weak field intensities ((w^) <IC {oJr)) the ion follows the flow as (vj^) ~ (v") 

• as the field gains in strength the ion starts to drift in a direction perpendicular to both the flow 
and the field until we get to the point where (v^) ~ when the field reaches high intensities 
{{ijjg) S> (wr))- The ion is then basically trapped in the field and its mean square velocity 
can be evaluated solely with equation (6), cr^ is found to be smaller than the flow velocity by 
a factor of a few. For example, if we choose Ai = 29 and An = 2.3 for the ion and neutral 
molecular mass numbers we find this factor to be ~ 3.5. That is, an ion would have on average 
an effective velocity at least ~ 3.5 times smaller than that of a neutral molecule of the same 
mass. 

This last observation has important implications. Let us assume that an observed emission line 
from a neutral molecular species is composed of contributions from a family of fiows (or eddies) of 
different velocities. We can then infer that for regions inhabited with a strong enough magnetic field 
which is, on average, not aligned with the local flow(s), we expect ionic lines to exhibit narrower 
profiles and a suppression of high velocity wings when compared to neutral lines. 



of the term cr?? and the fact that ( |v' 
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Obviously, we cannot expect all the flows to be perpendicular to the mean magnetic field. In 
any given case, there will be a distribution in the value of the angle existing between the directions 
of the flows and the field. We have only considered the extremities of this distribution (flows aligned 
or perpendicular to the mean field direction), in general for a given angle both velocity components 
(v|| and would have to be simultaneously included in the analysis. The total resulting effect 
observed will, in general, reside somewhere in between what is obtained for these special cases. 

It is also important to note that this phenomenon will only be observable in regions where the 
motions of the neutral particles do not have a zero mean velocity, i.e. (v") 7^ 0. For, if they did, 
equations (5) and (6) would have exactly the same form and the velocity dispersions would be the 
same for ions and neutrals. In other words, a thermal (or microturbulent) line profile would not 
show any manifestation of the presence of the magnetic field. 



2.2. Refined solution. 



It can be seen from equations (4), (6) and (9) that the amplitude of the effective velocity not 
only depends on the mean magnetic field strength but also on the amount of momentum transfered 
during a collision. It is therefore important to solve the problem for the more realistic cases where 
the collisions are not necessarily "head on". We, however, still retain the assumption that all are 
perfectly elastic. 
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The relaxation rate is now given by: 



(ur) = (acos(7))fc 



(15) 



a is the ratio of the change in velocity of the ion after a collision to the initial relative velocity 
between the two particles, 7 is the scattering angle of the ion as measured in its initial rest frame 
and the different averages ((a cos (7)), . . . ) are evaluated over the space of the scattering angle of 
the colliding neutral particle in the same frame of reference. We give in table 1 set of values for the 
different averages for two different ion molecular masses {Ai = 29 and 45). We would also like to 
point out that from numerical calculations we obtain the following relations (which could also be 
derived from simple physical considerations): 



(a cos (7)) ~ — 
rrii 

2 (a cos (7)) rrii 

which allow us to see that our first evaluation of the relaxation rate in equation (9) was wrong 
by a factor of two whereas the denominator in equation (14) for is unchanged. Another difference 
between equations (12)-(13) and (5)-(6) is that we now have a "mixing" between the two variances; 
the neutral dispersion parallel (perpendicular) to the magnetic field affects the ion velocity dispersion 
perpendicular (parallel) to the field. This is of course due to our more realistic treatment of the 
collisions where momentum can now be transfered between different directions. But despite these 
few changes, all the conclusions reached in the previous section still hold. 

In figure 1, we have plotted curves for the ion effective velocity (according to equations (10)- 
(15)) as a function of the mean magnetic field intensity when (^Vy*^ = 0, |(v^)| = 10 km/s and the 

neutral density n = 5 x 10^ cm~^. For simplicity, we also set (^Iv^p^ = (v")^. As can be seen, the 
field makes its presence felt even for relatively weak intensities {{B) > 10 fj,G) and the transition 
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between the regimes where ions follow the neutral flow to where they are trapped by the field is 
quite abrupt (1 /iG < {B) < 10 /xG). It is important to point out that it is, in principle, possible to 
evaluate the mean intensity of the magnetic field from this figure. One would have to compare the 
width of neutral and ion spectra and find the field strength that matches it at a given angle between 
the direction of the field and the line of sight to the observer. When this line of sight is parallel 

(perpendicular) to the field, the ion velocity should follow the curve for ^ (^^jj) ( (^^l))- There 

are, however, a few things that stop us from achieving this. Among these are the unknown amount 
of velocity dispersion inherent to the neutral flow and the uncertainty in the neutral density (used 
to determine the collision rate). We should also point out that recent observations indicate that 
fleld strengths in molecular clouds are in the range of a few hundreds ^G (Crutcher et al. 1999) 
where our curves in figure 1 are basically fiat and the ion velocity is insensitive to field strength 
variations. But we should not rule out the possibility that the determination of the field properties 
could become feasible in the future when using more complete models. Still at this point, it is 
nonetheless possible to assign a lower value for the field strength (given the density of the gas). 



3. Observational evidence. 

We would now like to bring support to our assertion of the previous section that we should 
expect that in certain conditions ion lines should have narrower profiles than neutral species. Fur- 
thermore, we contend that this situation is likely to happen frequently. The reason for this lies with 
equations (lO)-(ll). 

As we recall, the requirement is for a poor alignment between the local magnetic field and 
the "mean" flow; the effect being maximized when they are perpendicular to each other. Assuming 
that the magnetic field is strong enough, we can see from the aforementioned equations that ions in 
general will not follow the motion of the flow but will tend to drift away in a direction perpendicular 
to it. But since the total mean ionic velocity is minimal when the flow and the fleld are perpendicular 
to each other, charged particles will be more likely to aggregate in regions where this is the case. It 
is therefore very tempting to deflne two groups of objects (as far as the observation of the effects of 
the magnetic fleld is concerned): 

1. the field is aligned with the flow, no signiflcant differences between the spectra of comparable 
neutral and ion 

2. no general alignment between the field and the flow(s), the ion species exhibit narrower line 
profiles and suppression of the high velocity wings in their spectra. 

In what follows, we present observational evidence showing some examples that we believe belong 
to the second class. 



- 8 - 



3.1. Observations of HCN, HCO+ and N2H+. 

In HCN, HCO^ and N2H"'" we have three fairly similar molecules except for the obvious fact 
that the last two are ions whereas the first is a neutral. Indeed, they are all linear molecules with 
comparable atomic mass, similar rotation spectra and almost identical critical densities {ric — 10^ 
cm~^). They therefore appear to be excellent candidates to test our proposal. 

We present in figure 2 spectra of these three molecules in the J — > 3 — 2 transition for OMCl 
and J ^ 4-3 for OMC2-FIR4, OMC3-MMS6 and DR210H obtained at the Caltech Submillimeter 
Observatory with the 200-300 GHz and 300-400 GHz receivers. The observations (including the 
maps presented in figure 3) were done on several nights during the months of March, April, June 
and August 1999. Pointing was checked at regular interval using scans made on planets available at 
the time. Telescope efficiencies were calculated to be ^ 70 % for the 200-300 GHz receiver (beam 
width of ~ 32") and ~ 60 % for the 300-400 GHz receiver (beam width of ~ 20"). 

As can be seen from figure 2, all these spectra show the characteristics described earlier; the 
ion lines are all significantly narrower than the HCN lines and the suppression of the wings is also 
obvious. These conclusions are made more definite in table 2 where we present a comparison of the 
line widths (more precisely the standard deviations a,,) of the species for these molecular clouds. 
The widths were measured after the lines were modeled with a multi-Gaussian profile. We interpret 
the fact that the predictions made here on the differences in the appearance in the spectra between 
comparable neutral and ion molecular species seem to be easily observable in molecular clouds as 
strong evidence in favor of our assertions. 

We should however point out that the differences between the line profiles of neutral and 
ion species could be partly due to other factors. Indeed, the assumed coexistence of the different 
molecular species is not necessarily correct. Although we tried to choose ion and neutral species that 
are as similar as possible, it is likely that their observations sample different parts of the molecular 
clouds. In some cases, it might even be possible that two species are exposed to significantly 
different dynamical processes that could cause one of them to exhibit a narrower or larger line 
profile no matter what the effect of the magnetic field might be. Maps of the different molecular 
species for a given object could shed some light on this question. In the case of OMCl, Ungerechts 
et al. (1997) have extensively mapped the Orion molecular cloud (around OMCl) with 20 different 
chemical and isotopic molecular species. Amongst other things, they find an impressive degree of 
uniformity in the chemical abundances; but most interesting to us is how similar their maps of 
HCN and HCO"'" are ( J ^ 1 — 0, beam width of ~ 50"). On the other hand, N2H+ seems to have 
a very different spatial distribution (see their figure 2). We also present in figure 3 our maps of 
OMC2-FIR4 in HCN and HCO+ ( J ^ 4 - 3) done at the CSO. Although the ion species appears 
somewhat more extended, the two peaks are well aligned. This and the fact the two species have 
similar line center velocity (see figure 2) brings support to the coexistence assumption for these 
molecules in the 0MC2 cloud. 

The question of abundances is also important. For example, a given species could exhibit 
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more opaque lines which could change some characteristics of the line profiles (such as the relative 

importance of the high velocity wings, saturation or self-absorption). In fact, we believe that the 
differences between the spectra of N2H"'" and HCO^ are probably the result of this effect (see the 
case of DR210H in figure 2, the self-absorption feature at ~ —2.5 km/s is much weaker in the N2H"'" 
spectrum implying a lower abundance). 

Chemical differentiation could also play a role in explaining the differences between spectra. 
In this effect, we expect that our analysis only applies to long-lived ion species (like the protonated 
molecules considered here), other more reactive species (CO"^, SO"^, . . . ) should more or less behave 
as neutrals as almost every collision in which they are involved would entail a chemical reaction 
(Schilke 1999). 

4. Relaxation time and interaction between charged particles. 

In our analysis of the problem in section 2, we have considered the behavior of a given ion 
without taking into account the presence of any other charged particles. One might intuitively 
guess that for a weakly ionized plasma, like the ones probed with our observations at the CSO 
presented in the last section, it is probably safe to do so. 

One can make sure of this by comparing the mean time between collisions for charged particles 
in such a plasma with the time it takes an ion to relax to its steady state after being excited to 
a different velocity. The excitation could occur, for example, during a collision with an electron. 
However, because of the huge disparity between the masses of the two particles, it is unlikely that 
the ion would be much perturbed by such an encounter. A collision with another ion is more likely 
to produce a significant change in its velocity. 

At any rate, it turns out that the ion relaxation time is given by the reciprocal of the relaxation 
rate defined by equation (15), as one can verify when solving the equation of motion (1) for a mean 
collision force (Fc). The solution shows that the ion velocity decays exponentially with a time 
constant equal to (oj^)"^. When this time is compared to the mean collision time between ions, it is 
seen to be at least a few orders of magnitude smaller. Our simplified analysis is therefore justified. 

5. Conclusion. 

We argue that the presence of a magnetic field in a weakly ionized plasma can be easily detected 
through a comparison of ion and neutral line profiles. More precisely, we expect ion lines to often 
exhibit narrower profiles and significant suppression of high velocity wings. We have presented 
observational evidence obtained for four different molecular clouds which agrees with our theory. 

Because of the low intensities of field required for the effect to be noticeable, we expect the 
phenomenon to be widespread. However, where the mean magnetic field is aligned with the neutral 
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flow(s), or where the Hne width is dominated by the thermal width, should not show any significant 
differences between comparable neutral and ion spectra. These aspects of the problem will be 
treated in subsequent papers. 

We thank Prof. J. Zmuidzinas for his help during the preparation of this paper. M. Houde's 
work was done in part with the assistance of grants from FCAR and the Departement de Physique of 
the Universite de Montreal. The Caltech Submillimeter Observatory is funded by the NSF through 
contract AST 9615025. 



A. Derivations. 

Equations (3) and (4) can be easily derived by taking the mean of the equation of motion (1) 
and then assuming a steady state: 




Using equation (2) as an approximation for the force of interaction between ions and neutrals 
we get from the mean equation of motion: 



(Fcil) = (Al) 



e 



(Fc±) = -— (vx)x(B) (A2) 
ruiC 



with 



(F,) = -2A^e((v)-(v-)) (A3) 

where is the mean collision rate. From equation (Al) and (A3) it is straightforward to get 
equation (3) for (vy). Equation (A2) can be transformed to: 

(v±) = [((vx) - K)) X (CJ^)] 



into which we can insert equation (A2) using (A3) and finally get equation (4) which we rewrite 

here: 
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(o}^) and {ur) are defined in equation (8) and (9) respectively. 

Equations (5) and (6) for the velocity dispersions can be derived by using the following pro- 
cedure. We concentrate on the equation for a'j^. If we assume that the time of interaction during 
a collision is infinitively small, we can express v'^, the ion velocity (perpendicular to the field) 
immediately after a collision, as a function of the velocity just before the collision: 

= + ^ (A5) 
rrii 



with Ap± = — 2/x {v± — v" ) the change in the ion momentum. Upon taking the mean of the 
re of equati 
we are left with 



square of equation (A5) while once again imposing steady state conditions, i.e. ( Iv^^ ) = (|vx|^ 



AP±\ 
V. ) 



Apj 



(A6) 



Taking into account the fact that (v^) • (v] 
it is then straightforward to obtain: 



(vx) , as can be verified with equation (A4), 



rui _ -I 



(A7) 



The equation for ad follows from the same procedure. 



The more refined set of equations (10)-(15) can also be obtained in similar fashion but the task 
is rendered somewhat more complicated by the fact the interaction force Fc is not simply given 
by equation (2). We will omit their derivation here as they are somewhat lengthy and don't bring 
anything substantially new to the discussion. 

One last word concerning the energetics involved in the collision process. From the above 
discussion, it is apparent that an ion which starts off with the same velocity as the neutral flow 
would eventually settle into a steady state of lower kinetic energy (assuming a strong magnetic 
field). Once again focusing on the approximate head-on collision model as seen in the reference 
frame of the observer, we can understand this loss of energy by noting that the drag force felt by 
the ion (equation (2)) is proportional to the relative velocity between the colliding particles. As the 
ion is attempting to travel in a circular orbit around a given guiding center, this relative velocity 
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is greater when it is going upstream (and losing energy through collisions) than when it is going 
downstream (and gaining energy from the collisions). This will lead, on average, to a net loss of 
energy over a complete orbit. 
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Fig. 1. — Ion effective velocity (y ^v^y, y (v^) and a/ (v^)) as a function of the mean magnetic 
field strength when Vm' = 0, |v" | = 10 km/s, n = 5 x 10^ cm~^ and Ai = 29. 
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Fig. 2. — HCN (top), HCO"'" (middle) and N2H"^ (bottom) observations at the position of peak 
intensity of (clockwise starting from top left): OMCl, OMC2-FIR4, OMC3-MMS6 and DR210H. 
The bump at ~ 30 km/s in the middle spectrum of OMCl is a contamination from other species 
(SO2, i^CHaCN, ...). 
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Fig. 3.— HCN and HCO+ (J ^ 4-3) maps of OMC2-FIR4. Although the ion spatial distribution 
appears somewhat more extended, the two peaks are well aligned and the HCN and HCO"^ system- 
atic velocities agree (from figure 2). The two maps have the same lowest contour level of 5 K-km/s 
{'^ lOcj) and the following levels increase linearly with an interval of 5 K-km/s. The maps' grid 
spacing of 10" is half of the beam width (shown in the lower left corners). The pointing accuracy is 
better than ~ 5". 
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Parameter Ai = 29 Ai = 45 

(a cos (7)) 0.07734 0.05028 

(a2cos2(7)) 0.00778 0.00332 

(a2siii2(7)) 0.00359 0.00157 

(g^) 0.01137 0.00489 



Table 1: Numerical values for the different averages present in equations (10)-(15) for ion molecular 
masses of 29 and 45. 





Coordinates (1950) 


(T„(km/s) 


Object 


R.A. 


Decl. 


HCN 


HCO+ 


N2H+ 


OMC-1 


5 32 47.2 


-05 24 25.3 


17.42 


3.23 


1.87 


OMC-2 FIR 4 


5 32 59.0 


-05 11 54.0 


3.59 


2.72 


1.34 


OMC-3 MMS 6 


5 32 55.6 


-05 03 25.0 


1.40 


0.71 


0.56 


DR 21(0H) 


20 37 13.0 


42 12 00.0 


5.75 


4.61 


2.08 



Table 2: Comparisons of line widths (standard deviations a^) between HCN, HCO''", and N2H"'" for 
the four objects presented in figure 2. 



